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Abstract 

We study the formation of electron-hole pairs for disordered systems in the 
limit of weak electron-hole interactions. We find that both attractive and 
repulsive interactions lead to electron-hole pair states with large localization 
length A2 even when we are in this non-excitonic limit. Using a numerical 
decimation method to calculate the decay of the Green function along the 
diagonal of finite samples, we investigate the dependence of \2{U) on disorder, 
interaction strength U and system size. Infinite sample size estimates £,2{U) 
are obtained by finite-size scaling. The results show a great similarity to 
the problem of two interacting electrons in the same random one-dimensional 
potential. 
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It is well-known that randomness leads to localization of non-interacting electrons and 
holes. This effect is especially strong in one-dimensional (ID) systems leading to com- 
plete localization even for very small disorder ||^. On the other hand, in the limit of strong 
attractive Coulomb interaction, electrons and holes will pair into excitons with large life- 
times. In the present work, we concentrate on the intermediate problem of weakly interacting 
electrons and holes (lEH) in a random environment. The energy scales are such that the 
band width is larger than the disorder W which in turn is larger than the interaction strength 
U and so we do not have bound excitons, but rather electron-hole pairs. Such a problem 
is relevant for the proposed experimental verification of the two-interacting particle (TIP) 
effect by optical experiments in semiconductors . The TIP problem has recently attracted 
a lot of attention after Shepelyansky argued that attractive as well as repulsive on- 
site interactions between two bosons or fermions in a single random potential lead to the 
formation of particle pairs whose localization length A2 oc U'^Xf is much larger than the 
single-particle (SP) localization length Ai. His major prediction is that in the limit of weak 
disorder a pair of particles will travel much further than a SP even for repulsive interaction 
U >0. 

Here, we consider the effect of onsite interaction on a single electron-hole pair, modelled 
by TIP in two different ID random potentials. The Hamiltonian we consider is 

H = —t (|n, m){n -|- 1, m| -|- \n, m) {n, m + 1\ + h.c.) 

n,m 

+ ^\n,m) {e\ + el^ + U5nm) {n,m\ (1) 

n,m 

where for the case of TIP in different potentials, e.g. two electrons on neighboring chains, 
or an electron and a hole on the same chain, we have e\ 7^ e^. For simplicity, both e\ 
and are chosen randomly from the interval [—W/2,W/2\. In the following we call this 
situation the lEH case. We will show that the results for lEH are similar to the (standard) 
TIP problem when both particles are in the same potential e\ = e\. In both cases U is the 
on-site interaction between the two particles. We use hard-wall boundary conditions and the 
hopping element t = 1 sets the energy scale. We note that this Hamiltonian corresponds to 
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the 2D Anderson model if we replace e^ + e^ in ([l|) by e„m and choose e„m £ [—W/2, W/2]. In 
this case U is an added on-site potential. This has been recently used to test our numerical 
method by comparing with more established methods valid for this model |^. 

To obtain our results we use a decimation method (DM) This involves replacing 

the full Hamiltonian by an effective Hamiltonian for the doubly-occupied sites only. It 
should be stressed that this method is exact and no approximations have to be made in the 
decimation process. It is then possible, via a simple inversion, to obtain the Green function 
matrix elements (1, l|G2|Af, M) between doubly-occupied sites. We shall be focusing on the 
lEH localization length A2 obtained from the decay of the transmission probability of lEH 
from one end of the system to the other. This is defined by 

;^ = -p^^ln|(l,l|G2|M,M)|. (2) 

In order to reduce possible boundary effects, we compute A2 by considering the decay between 
sites slightly inside the sample. We present results for the band centre, i.e., energy E = 
for 14 disorder values W between 0.5 and 7, for 21 system sizes M between 51 and 251, and 
11 interactions strengths U = 0, 0.1, . . . , 1.0. For each triplet of parameters {W, M, U) we 
average the localization lengths A2 computed from the Green function according to Eq. 
over 100 samples. 

In Fig. |l], we show the lEH results for M = 201. Let us first turn our attention to the 
case f/ = 0. As pointed out previously |jTO|, the TIP Green function G2 at = is given 



by a convolution of two SP Green functions Gi at energies Ei and —Ei. This implies that 
A2 ~ Ai/2. We have therefore included data for Ai/2 in Fig. [l| which we have computed 
by a transfer matrix method (TMM) |11]] in ID with 0.1% accuracy. Comparing these 
results to the localization lengths A2 obtained from the DM, we find that for 1 < < 6, 
the agreement between X2{U = 0) and Ai/2 is rather good and, contrary to TMM results 



there is no large artificial enhancement at f/ = 0. For smaller disorders W < 1, 



we have A2 ~ M/2 so that it is not surprising that the Green function becomes altered due 



to the finiteness of the chains [0. This results in reduced values of A2. It is noticeable 



from these results, however, that the values of X2{U = 0) are still slightly larger than Ai/2. 



This is similar to TIP results |T0| and to a numerical convolution of the SP Green functions 
calculated by exact diagonalisation [0. For U between 0.1 and 1 and W > 1.2 we have 
found that the localization lengths X2{U) are increased by the onsite interaction as shown in 
Fig. |I]. For [/ ^ 1, it is well-known that the interaction will split the single TIP band into 
upper and lower Hubbard bands and, hence, we expect that for large U the enhancement of 
the localization length at E = will vanish. In Fig. |^ we present data for X2{U)/X2{0) for 
U = —4, . . . , 4. We first observe that at the band centre the enhancement is symmetric in 
U. This is why we usually only consider f/ > in agreement with the previous arguments 
and calculations for the TIP case [|7|,P, p!0| , p!2| -p!7[] . We have checked that away from the band 



centre the enhancement is asymmetric in U. For small \U\, we see from Fig. Q that the 
localization length increases nearly linearly in |f/| with a slope that is larger for smaller W 
and we do not see any f/^ behavior as argued in Refs. PJ6|JT8||. At large \U\ the enhancement 
starts to decrease again. For TIP, it has been suggested that there exists a duality for 
U and v^24/f/ for very large \U\. The crossover between the two asymptotic regimes should 
accor at Uc = 24^^^. For lEH, we find that within the accuracy of our data, we can argue for 
an agreement with the duality. As for our TIP data |^ , we observe the best lEH agreement 
with duality for W = 5 but the maximum enhancement maxjj [X2{U) / X2{0)] still seems to 
depend upon the disorder. 

In order to overcome the problems with the finite chain lengths, we construct finite-size 
scaling (FSS) curves for each U and compute from these scaling parameters which are the 
infinite-sample localization lengths C,2{U). This method has been proven very useful for the 
non- interacting case |]19| and recently for TIP studies ||7|,p!0|. In Fig. ^ we show the raw data 



of the reduced lEH localization lengths X2/M which are to be scaled just as in the standard 



TMM |T9|. Note that data for small W are rather noisy and will thus most likely not give 
very accurate scaling. In order to set an absolute scale in the FSS procedure, one usually fits 
the smallest localization lengths of the largest systems to X2/M = ^2/M + 6(^2/^)^ with b 



small [0. Due to numerical problems of estimating a small localization length of the order 



of 1 in a large system by Eq. (^) we instead fit for each U to the locahzation length atW = 3 
and adjust the absolute scale of ^2 accordingly. In Fig. ^ we show the resulting scaling curves 
A2/M = f{C,2/M) for U = 0, 0.2 and 1. The above mentioned numerical errors in the data 
at large M and W are visible only in very small upward deviations from the expected 1/M 
behavior. The results are very similar to the TIP problem. It is interesting to note that it 
is even possible to scale the present lEH data together with the data previously obtained 
for the TIP case 0. From this more accurate scaling we compute the scaling parameters 
^2 which we show in Fig. ^. The power-law fits ^2 oc W'"^^ to the data with W G [1,5] 
yield an exponent a which increases with increasing U as shown in the inset of Fig. e.g., 
a ^ 1.1 for [/ = and a ^ 1.5 for U = 1. Thus, although in Fig. |I] the A2 data at M = 201 
nicely follows Ai/2 for U = 0, we nevertheless find that after FSS with data from all system 
sizes, ^2(0) as against Ai/2 still gives a slight enhancement. Because of this we will in the 
following compare ^2{U > 0) with ^2(0) when trying to identify an enhancement of the 
localization lengths due to interaction. For comparison, the exponents obtained from the 
same fit applied to the TIP problem [[^ are also shown. Note that, as expected from Fig. |^, 
FSS is not very accurate for small W. Therefore, in what follows we shall only use ^2 values 
obtained ior W > 1. For data corresponding to < 1, we actually used the extrapolated 
values of ^2 from the power-law fit to continue the FSS curves of Fig. ^ to < 1. 

We now compare our lEH results with various fits proposed for TIP. From an effective 
random matrix model P JT5| , p!6[| A2 oc A^ was obtained for large values of Ai. To correct for 



smaller values of Ai a more accurate expression was suggested [jT6[ to be A2 oc Af (1 + c/Ai). 
It is important that /3 in this work depends on U and ranges from 1 at small U and very 
large U to nearly 2 ior U ^ t. As discussed above we translate this fit function into 
C,2{U) oc ^2(0)'' ^1 + III Fig- H we show respective data for disorders W G [1,6]. The 

fits are good and reflect in particular the deviations from a simple power-law ^2iU) oc ^(0)^ 
for small localization lengths. In the inset we present the dependence of (3 on U: we find 



P < 1.5 for all U values considered unlike Ref. |[T^. The values obtained for the TIP case 
are shown for comparison. 



In Ref. [H the functional dependence of the TIP localization lengths A2 = Ai/2 + c|f/|Af 
has been suggested. Taking instead of Ai/2 the more suitable ^2(0) this can be translated as 
6(f/) -6(0) oc \UMO?. In Fig. 0we plot [6(f/) - 6(0)]/|f/| vs. g{UMO) for U G [0.1, 1] 
where we have chosen g{U) so that the data for different U can be placed on top of the 
U = 0.1 data. In the inset of Fig. |^ we show that g{U) starts to deviate from 1 already for 
U > 0.4. Thus we see that the linear behavior in \U\ found in Ref. [§] holds only for very 
small U in agreement with TIP. We obtain a good fit to [^(f/) — 6(0)]/|t/| with a single 
exponent /5 = 1.6 ± 0.1 instead of 2. This is somewhat different from the fits for TIP which 
give /5 ^ 2 for small 6(0) < 10 and /3 ~ 3/2 for larger 6(0) 0- "^^^ reduction of the slope 
below j3 = 2 may be due to insufficient disorder averaging and thus an underestimation of 
the FSS results. 

In conclusion, we have presented detailed results for the localization lengths of electron- 
hole pair states which may be realized in the non-excitonic limit of optically excited semicon- 
ductor heterostructures |^. We observe an increase of the two-particle localization length 
due to onsite interaction in the band centre. This suggests the formation of an electron-hole 
pair with possibly enhanced transport properties. We emphasize that our results apply to 
the non-excitonic limit with bandwidth larger than W > U. We have fitted our data to var- 
ious suggested models with varying success. The results are all similar to the standard TIP 
problem in a single random potential and thus we conclude that the case of an interacting 
electron-hole pair is very close to the TIP problem |^. 
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R.A.R. gratefully acknowledges support by the Deutsche Forschungsgemeinschaft through 
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FIG. 2. Enhancement A2(?7)/A2(0) as a function of interaction strength U at E = iov 
disorder W = 3 (+), 4 (*), and 5 (x) and M = 201. The thick (thin) hnes indicate data for 
U > {U < 0), fuh (open) circles denote the maximum for each disorder. The dashed hne marks 
Uc = 24^4. 
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FIG. 3. Reduced localization lengths X2/M for U = (left) and U = 1 (right) for 14 disorders 
between 0.5 (A) and 7 (N) as in Fig. |. 
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FIG. 4. Finite-size scaling plot of the reduced localization lengths A2/M for various U. The 
data for U = 0.2 (1) have been divided by 2 (4) for clarity. Data corresponding to W = 1 are 
indicated by filled symbols. 
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FIG. 5. Localization lengths ^2 after FSS for various U. The solid line represents Ai/2, the 
dashed lines indicate power-law fits. Inset: Exponent a obtained by fitting ^2 W~'^°^ to the 
data for each U (filled symbols). The result for the standard TIP problem is also shown (open 
triangles) . 
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FIG. 6. (,2(U) after FSS for various U plotted versus ^2(0). The lines are fits 
^2{U) oc ^2(0)^[1 + c/^2(0)]. Inset: Exponent P from the fits (filled symbols) and for the stan- 
dard TIP problem (open triangles). 
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FIG. 7. Enhancement of the locahzation lengths [£,2{U) — S^2{^)]/\U\ plotted for W G [1,7] ver- 
sus 5(^7)^2(0) where g{U) was obtained by a mean-least-squares fit procedure to make all the data 
compatible with g{0.1) = 1. The straight line is the curve ^(f/) = 6(0) + 0A9\U\[g{U)^2{0)V-^^ ■ 
Inset: Behaviour of g{U) for lEH (filled symbols) and for TIP (open triangles). 
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